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NOTATION 

A  Source  strength 

b  Exponent  of  density  stratification 

f  Depth  of  submergence  of  disturbances 

g  Acceleration  due  to  gravity 

H  Depth  of  thermocline 

M  Quadrupole  strength 

N  Vaisala  frequency 

t  Time 

u,v  Velocity  components  in  x,y  direction  respectively 

U,V  Transformed  velocity  components 

i 

x,y  Cartesian  coordinates 

a  Horizontal  component  of  wave  number 

| 

U)  Frequency 

p-  Equilibrium  density 

t  Thickness  of  thermocline 

6(x)  Dirac  delta  function 

cp  Fourier  transform  of  V 
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INTRODUCTION 

A  submerged  disturbance  in  a  stratified  ocean  will,  in  gen¬ 
eral,  generate  motions  that  are  quite  different  from  those  gen¬ 
erated  in  a  homogeneous  ocean.  By  stratification  we  mean  a  den¬ 
sity  stratification  due  to  either  the  presence  of  a  salinity  or 
temperature  gradient  or  both,  the  fluid  being  assumed  to  be  in¬ 
compressible.  The  difference  is  not  simply  a  modification  of 
the  motion  by  the  density  stratification,  which  is  slight  in  a 
real  ocean,  but  rather  the  appearance  of  other  forms  of  motion 
whose  existence  depends  entirely  on  the  density  stratification. 
The  difference,  then,  is  one  of  kind  and  not  simply  one  of 
degree . 

In  this  report,  the  surface  currents,  i.e.,  currents  at  the 
air-sea  Interface  due  to  a  two-dimensional  submerged  disturbance, 
are  considered.  The  equilibrium  density  profile  of  the  ocean  is 
assumed  to  be  of  the  form  shown  in  Figure  la.  This  profile  is 
chosen  for  its  simplicity  and  also  for  its  close  approximation 
to  reality.  The  disturbance  is  assumed  to  be  always  below  the 
thermocline:  It  is  a  quadrupole  of  constant  strength  amplitude 

and  with  a  simple  harmonic  time  dependency.  In  spite  of  its  sim¬ 
plicity  it  provides  a  time  scale  (inverse  frequency)  and  length 
scale  (proportional  to  its  depth  below  the  thermocline)  charac¬ 
terizing  the  disturbance. 

In  the  Appendix  similar  calculations  are  carried  out  for  an¬ 
other  class  of  density  profiles  shown  schematically  in  Figure  lb. 
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FORMULATION  OF  THE  PROBLEM 


The  equilibrium  density  profile  of  the  ocean  is  taken  to 
be  of  the  following  form: 

/  , 


?-< 


Pi 


Pa 


-2b!  (y-f ) 


in  R  :  (f  +  H>y>f  +  T) 


in  Rx  :  (f  +  T  >  y  >  f ) 


e  2t>3y  in  Rs  :  (f  >  y  >  -  ®) 

The  -fluid -is -assumed  to  be  inviscid  and  incompressible.  In  the 
above,  p,  px ,  pa ,  f,  t  and  H  are  all  positive  real  constants.  A 
schematic  diagram  of  the  profile  is  shown  in  Figure  la. 


We  are  interested  in  the  currents  at  the  free  surface 
y  =  f  +  H  generated  by  a  two-dimensional  quadrupole  of  constant 
strength  amplitude  M  and  with  a  simple  harmonic  time  dependency 
of  frequency  ou  located  at  the  origin  of  the  coordinate  system 
(x  =  0,  y  =  0).  The  solution  for  a  periodic  source  of  strength 
A  is  found  first.  From  this  solution,  the  one  for  a  quadrupole 
is  obtained  by  differentiation. 

It  is  convenient  to  work  with  the  transformed  velocity 
components  U  and  V  defined  by 


U  =  p  ?u 
V  =  p  ^v 
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where 

u,  v  are  the  physical  velocity  components. 

Since  the  forcing  function  is  simple  harmonic,  so,  after  a 
sufficient  time  has  elapsed  for  the  transients  to  die  away,  the 
whole  system  would  acquire  a  simple  harmonic  motion  of  the  same 
frequency.  In  this  final  quasi -steady  state,  the  time  dependency 
of  the  dependent  variables  may  be  factored  out  so  that 


U(x,y, t)  =  U(x,y)e~ 
V(x,y,t)  =  7(x,y)e"iU)t 


If  the  disturbance  is  weak, so  that  the  resulting  motion 
may  be  considered  as  small  perturbations  about  the  state  of 
hydrostatic  equilibrium,  then  the  differential  equations  gov¬ 
erning  V  in  the  different  regions  are  (see  Wong,  1965) 


d2  V  d2  V 

in  R  - -  +  - -  =o  [  l] 

0  By2  dx2 


in  Rx 


by3  1  uu2  I  dx2 


bi2  Vx 


0 


in  Ra 


+ 1 1  -  -  ba2 v3 

dy2  '  W2  /  dx2 


[2] 


2i r /"p  A  6 (x)  [b(y)  -  ba6(y)] 


[31 
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where  B  denotes  the  Dirac  delta  function  and  a  dot  denotes  a 
differentiation.  Subscripts  have  been  added  to  V  to  denote  the 
function  in  the  various  regions.  The  right  hand  side  of  Equa¬ 
tion  [3]  represents  the  source. 

The  surface  current  is  essentially  due  to  the  formation 
and  propagation  of  internal  waves.  Because  of  the  much  larger 
density  discontinuity  at  the  air-sea  interface  the  vertical 
movements  there  are  negligible  for  all  practical  purposes. 
Therefore,  the  boundary  condition  there  can  be  taken  to  be 


7=0  at  y  =  f  +  H 
o 


[4] 


At  y  =  f  and  y  =  f  +  T  both  U  and  V  have  to  be  continuous. 
This  gives  rise  to  the  following  boundary  conditions: 


V  =  Vj.  at  y  =  f  +  t 
o 


7X  =  Va  at  y  =  f 


[5] 

[6] 


&V 

+  at  y  =  f  +  T 


[7] 


at  y  =  f 


[8] 
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The  last  two  equations  were  obtained  from  the  equation  of  con¬ 
tinuity,  viz., 

|f  +!f+bV=2Tr/pA  6(x)  6(y) 

which  becomes 


SUo  3Vq 
dx  +  dy 


0 


in  R 

o 


[93 


0 


in  Rx 


ClO] 


3x  + 


SV2 

3y 


+  b2  V2 


2tt  /p  A  6  (x)  6  (y ) 


in  Rg 


cii] 


In  addition  to  the  above  conditions  other  conditions  have 
to  be  imposed  at  Infinity  to  insure  a  unique  solution.  These 
additional  conditions  will  be  discussed  as  they  arise. 


SURFACE  CURRENTS  DUE  TO  A  SUBMERGED  SOURCE 


Let  cp^  denote  the  Fourier  transform  of  with  respect  to 
x  so  that 


i 


+  » 


i  =  0,  1,  2. 
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Then,  from  the  equations  governing  the  following  set  of 
tlons  for  Is  obtained: 

cpQ  -  a3cpQ  =  0  for  f  +  H>y>f  +  T 


cpi  - 


bn3  + 


L  _  2gb,J 

Q3 

1  U)3  J 

cpx  =  0  for  f  +  T  >  y  >  f 


cpa  - 

b  2  +(x  _  2gbs_1 

jo? 

1 

l  3  \  u>3  J 

cp2  =  2v  A  /p[6(y)  -  b6  (y  )J  for  y 


cpQ  =  0  at  y  =  f  +  H 

cpQ  =  cpi  at  y  =  f  +  T 

cpi  =  cpa  at  y  «  f 

cpQ  =  bicpi  +  cp!  at  y  =  f  +  T 

bx q>!  +  cpi  =  bacpa  +  cpa  at  y  =  f 


equa- 


[12] 

[131 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 


The  most  general  solutions  of  Equations  [12]  to  [ 14] 
are  respectively. 
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-  a|y  + |a|  y 


cp  =  A  e  +  B  e 

oo  o 


[20] 


-Pi  y  +0i  y 

cpx  =  Ai  e  +  B].e 


[21] 


-03y  +02  y  ih  »  -03y 

cpg  =  A3e  +  Bge  +  H(y)  w  A  /pj^-  +  lj  e 


+  [1  -  H(y)].ir  A  /p 


* 


+03  y 


[22] 


where  A^,  are  arbitrary  functions  of  the  parameter  a,  H(y) 

denotes  the  Heaviside  function  and  fa ,  fa  are  given  by 


Because  of  Equation  [15]>  cpQ  must  be  of  the  form 

cpQ  =  Aq  sinh  | a)  [(f  +  H)  -  y]  [23] 

Aq  now  is  some  other  arbitrary  function  of  a.  Because  of  condi¬ 
tions  at  y  -*  -  Aa  s  0  and  Equation  [22]  becomes 
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cPs 


Ebe+^ay  +  H(y)  tt  A  /p  (  ^  +  1 

Pa 


e-^y 


+  [  1  -  H(y )  ]  tt  A  /p 


biL 

03 


+0a  y 


[24] 


Substituting  known  expressions  for  cp^  into  Equations  [ l6] 
to  [19]  yield  the  following  set  of  four  equations  for  the  deter¬ 
mination  of  Aq,  Ax  ,  Bi  and  Bz  : 


Aq  sinh  | a |  (H  -  T)  -  Axe-*<f+T>  +  Bxe+^f+T> 


|a|A  cosh  |a|  (H  -  t)  =  Ax  (bx  -  0x )e“^x (f+T Bx  (bx  +px  )e^x  ^f+T  ^ 


Ba  e^s  f  +  tt  A  /p  j 


0a 


+  1 


-03  f  A  -0i  f  _  +fi,  f 
e  =  Ax  e  r  +  Bx  e  r 


Eb  (ba  +  )e^3  f+  (b3  -  02  )  ir  A  /"p 


j£-  +  lle^2f  =  A,  (bx-Px  )e‘^f 


+  Bx  (bx  +Px  )e+^x  ** 


Solving  the  above  equations  simultaneously, Aq  is  found 
to  be : 

A _ 4ttA  /p  (ba+  0a  )e~^af _ 

°  (?+R)(bx-ba-0i-p3)e+^lT+  (?  -n)  (bx  -ba  +0i  -0a  )e^1  T 


[253 
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?  =  sinh  jet f  (H  -  t) 

ri  =  [  |  a|  cosh  |a|  (H  -  t)  +  bj.  sinh  )  a  I  (} 

Pi  1 


Taking  the  inverse  transform,  Vo  is  expressed  as  a 


integral 


.+» 


V  =  —■  I  A  slnhlal  [ (f  +  H)  -  y]eiaX  da 
o  2ir  o  ii' 


U  can  be  found  from  Equation  [9]  as  follows: 

0 


dU  dV 

_ o  _  _ o 

dx  “  ”  by 


.  +  » 


=  + 


2ir 


A  |a  |  coshjaj  [  (f  +  H)  -  y]e 


iax 


U 

0 


C(f  +  E)  -  ylalM 


In  arriving  at  the  last  expression  for  UQ  the  fact  that 
antisymmetric  in  x  was  used. 


:  -  t)] 

Fourier 


[26] 


da 


da  [271 


U  is 

O 
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The  surface  current  Induced  by  the  source  is,  then, 
given  by 


u 

o 


-iuit 

e 

2TTl/pi 


[28] 


In  what  follows  we  shall  always  assume  that 


e  = 


£2  -  Pi 

p3 


«  1 


First  consider  the  simplest  special  case  where  the  thermo- 
cline  is  a  discontinuity  (t  =  0)  and  the  layer  below  the  thermo - 
cline  is  homogeneous  (b2  =0,  p  =  pa ) „  For  this  special  case 


A  = 
o 


2 tt  A  /p^  a  e 


-  |  °t  |  f 


+  a  H  a 

|  a  |  e  -  sinh  |a|  H 


U)2 


Therefore 


u  =  -  i  A  e 
o 


-iuut 


I-*'  ^Me~|aKt+H)elaXda 

-L os  ou2  -  ge|a|e~la^  sinh|a|H 


[29l 
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Since  the  integrand  has  poles  on  the  real  axis,  the  Integral  is 
undefined.  Under  such  circumstances  it  is  usual  to  deform  the 
path  of  integration  so  that  it  goes  around  instead  of  through 
the  poles.  How  the  path  is  to  be  deformed  depends  on  additional 
physical  information.  Here,  the  surface  current  is  a  conse¬ 
quence  of  the  formation  and  propagation  of  interfacial  waves  at 
the  thermocline.  These  waves  are  all  outgoing  waves  from  the 
central  region  of  disturbances.  Therefore,  the  current  also 
possesses  this  outgoing  wave  property.  For  such  a  situation  the 
path  has  to  be  deformed  slightly  above  the  real  axis  when  x  <  0 
and  slightly  below  when  x  >  0.  The  path  of  integration  can  be 
closed  by  a  large  semicircle  in  the  upper  half  plane  when  x  >  0 
and  the  lower  half  plane  when  x  <  0.  It  can  easily  be  shown 
that  the  contribution  from  the  Integration  along  the  semicircle 
vanishes  as  the  radius  tends  to  infinity.  Since  the  integrand 
is  analytic  except  for  poles, the  integral  may  be  evaluated  by 
the  method  of  residues.  The  contributions  from  the  poles  having 

a  non-vanishing  imaginary  part  will  vanish  exponentially  with 

ictx 

x  because  of  the  factor  e  .  Therefore,  for  x  large  and  posi¬ 
tive  only  the  positive  real  roots  of  the  equation 

uua  -  gea  e~  sinh  aH  =  0  [30] 

are  important.  This  equation  is  plotted  in  Figure  2.  It  is 
seen  that  for  a  given  ui  there  is  only  one  real  positive  a, 
say,  aQ  which  satisfies  the  above  equation. 
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i 


Hence,  for  large  and  positive  x 

-a  (f+H)  . , 

,  o  i(a 


U  ~  -  7T  A 

o 


u)  e 


-  U)t ) 


do 

da  , 

\  la 


where 


"lalH 

a8  =  ge  |a  |e  sinh  |a  |H 


C  3 1] 


and  aQ  is  the  positive  real  solution  of  Equation  [30]. 

Next  consider  the  effects  of  a  thin  thermocline  on  the  sur¬ 
face  currents  generated  by  the  source,  i.e.,  T  is  small  compared 
with  H  but  finite.  As  in  the  previous  case  b3  is  set  equal  to 
zero.  Now 


i 


i 

i 

i 


bx  =  ^  m  ( 1  +  « )  w 


therefore 


Pi  T  -  T-Wbi8  + 

[l  2sbl 

1  * 

V 

i  U)2 

1 

fe3  ~~  3 

o 

.  a3 

=v~ 

t3  -  g  e 

uu3 
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If  t  Is  small,  its  effects  on  the  dispersion  relation  (Equa¬ 
tion  [30])  may  be  expected  to  be  small.  The  third  term  under 
the  square  root  can  then  be  approximated  by 


t  a _ 

-aH  .  , 
e  slnh  aH 


1 


2ta 

-2aH 

e 


f  2Ta 


- 


for  aH  large 


t 

—  for  aH  small 
H 


From  Equation  [31]  it  is  seen  that  only  waves  with  small  wave 

numbers  produce  significant  currents  because  of  the  factor 
- 1  a  I ( f +H ) 

e  .  Therefore,aT  is  small  since  t/H  Is  assumed  to  be 

small.  Hence,  /3X  t  is  small  and  we  have 


e^1  T  «  1  +  j3x  T 
e~0l  T  «  1  -  /3i  T 

Under  these  circumstances  and  with  ba  =  0,  Aq  simplifies  to 


_  -|a|f 

_ 2 TT  A  /Tte  la  je  _ 

0  |  ot  |  ( 1  -t-  |a|  T  )e~  ia  lH  -  ge  —  sinhj  a|  H 
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and  the  surface  current  to 


+  ® 


u  =  -  1  A  e 

o 


-lout 


J  It.  I  -I“l <f+H-T) 


g(l  4-  a  T ) 


lax  , 
e  da 


.  -lal(H-T) 

-  ig+|a|T  e  slnh | a | (H-T ) 


C32] 


The  Integral  In  Equation  [32]  can  be  evaluated  as  In  the  previous 
case  where  the  thermocline  Is  a  discontinuity.  For  large  posi¬ 
tive  x  the  surface  current  is  given  by 


u 


o 


-ttA 


-ao(f+H-T)  l(aox-uut) 


(l  +  a  t]  I 

(dyl 

l  0  j 

[da  1 

a 

o 


[33] 


where 


_SL| 

a 

1+1 

a  T 

-  |a  1  (H-t  ) 

e  sinhlal(H-T) 


and  aQ  is  the  positive  real  root  of  the  following  equation: 


gea 

1  +  a  t 


-a(H-T) 

e  sinh  a(H-T) 


=  0 


[34] 


Finally,  let  us  consider  the  effects  of  a  stratification 
below  the  thermocline;'' the  thermocline  is  here  taken  to  be  a 
discontinuity,  i.e.,  T  =  0,  ba  is  small  so  that 
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&  -  V1*3  +I1  -  ^Ha2  *  M  V1  "  ^ 


Under  such  circumstances  Aq  simplifies  to 


'irA/p  "\/l  - 


2%e 

UT 


i°i r 


1  -  sinh  !  ct  |  H  +  cosh  |a|  H  -  -  -  sinh  |  at  |  H 


and  the  surface  current  to 


u  «  -  i  A  e 

o 


-itl)t  r - -r-r- 

L  2gb 

1~dri 


M 


-V1  l“lf 


a 


iax 

e  da 


sinh  |  a|  H  +  cosh  |  a|H  -  S^Ja'  slnh|a|H 


[351 


When  t U3  <  2gb2 ,  the  Integrand  has  no  pole  on  the  real  axis,  This 
can  be  shown  as  follows:  Let  aQ  be  a  real  root  of  the  equation 


!  1 


:"2gbL 


sinh  |a|  H  +  cosh  |a|  H  -  sinh  |aj  H  =  0 


[36] 
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If  it ?  <  2gba  then  aQ  has  to  satisfy  the  following  pair  of 
equations: 

~  1  slnh  |o0|  H  =  0 


cosh  la  I  H 

1  o 1 


slnh  |aQ|  H  =  0 


but  this  is  impossible.  Therefore,  Equation  [36]  has  no  real 
solution  for  uua<  2gba  .  However,  for  every  uf  >  2gba  there  is  a 
unique  positive  solution  (see  Figure  2).  Let  this  be  denoted 
by  aQ.  Then  for  large  positive  x  the  surface  current  is  given 
by 


2ttA 


2gba 

uu3 


a 

o 


i  (a  x  -  out ) 
e 


[37] 


where 

C  =  ~\Jl  -  sinh  |a|  H  +  cosh  |a|  H  -  sinh  |a|  H 


SURFACE  CURRENTS  DUE  TO  A  SUBMERGED  QUADRUPOLE 

The  quadrupole  under  consideration  has  a  configuration  shown 
in  Figure  3 .  The  strength  M  of  the  quadrupole  is  given  by 

M  =  A3  A3 
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where  A  is  the  strength  of  the  individual  sources  and  sinks  and 
A  is  the  distance  separating  them. 

The  surface  current  w  due  to  a  quadrupole  of  strength 
amplitude  M  and  time  dependence  located  at  the  origin  is 

obtained  by  differentiating  the  source  solution,  i.e.,  u^  with 
respect  to  x  and  replacing  A  by  M. 

Let  us  choose  H  as  the  reference  length  and  /H/ge  as  the  ref¬ 
erence  time  and  introduce  the  following  dimensionless  variables 


t'  =  t/Vll/ge 
x'  =  x/H 
f'  =  f/H 
T'  =  T/H 

a  '  =  a  H 

o  o 


The  surface  currents  generated  by  the  quadrupole  corre¬ 
sponding  to  the  three  cases  analyzed  in  the  previous  section 
for  the  source  are  found  to  be 
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Case  I  t  1  =  0,  b3  =  0 

For  large  positive  x  the  surface  current  due  to  the 
quadrupole  Is  given  by 


W 


ttM 

w 


ui'(a  ')8  e  ° 


a  '(f'+l)  l(a  'x'-ou't') 


i  da '  \ 
’  da' 


a 


The  magnitude  of  the  current  is 


where 


M 


ttM 

iF 


-  a  ' (f '  +  1) 
iul(ao1)8  e  0 _ 

/da'] 

da'  . 
i  /  a  ' 
o 


U* 


(a')a  =  |  a'  |  e  sinh  |  a'  | 


[38] 


and  aQ'>  which  is  the  asymptotic  wave  number  of  the  current,  is 
the  positive  real  root  of  the  equation 


(u>')a  -  a'  e” 


sinh  a' 
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Case  Ii  T 1  4-  0,  b2  =  0 

For  large  positive  x  the  surface  current  due  to  the 
quad rup ole  is 


.  .  id  '  ( a  '  )s  e  0 

1  irM  v  o  ' 

w~  iF 


-a  '(f'+l-T')  i (a  'x'-uj't'  ) 


(1  +  a  « T  • ) 
x  o  ’  da1 


a 


The  magnitude  of  the  current  is 


I  w  |  ~  |  jp- 


-a  '(f'+l-T') 


(1  +  VT,)  (to'j  , 

l  /  a  ' 


irM 

iF 
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C393 


where 


-|a'|(l-T') 

(t  '  )a  =  -jy  |  T-Ty  e  sinh  I a '  I  ( 1  -  t  ' ) 


and  a  1  is  the  positive  root  of 


ID1 


a1 


-a '  ( 1-T  '  ) 


1  +  a'T' 

When  T '  *  0,  this  case  reduces  to  Case  I. 


sinh  a' (1  -  T' )  =  0 
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Case  III  t'=  0,  b3  ^  0 


For  large  positive  x  the  surface  current  due  to  the 
quadrupole  is 


,  ttM 

w  ~  r 
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u* 


where 


C  1  =  yl  —  ( ~r)  sinh|a’  |  +  cosh|a'  |  -  sinh|a' 


and  aQ 1  is  the  positive  real  root  of  the  equation 


C  =  0 


as  N'  tends  to  zero  this  case  tends  to  Case  I. 
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DISCUSSION  OF  RESULTS 

The  surface  current  due  to  a  two-dimensional  quadrupole 
of  constant  strength  amplitude  and  simple  harmonic  time  de¬ 
pendency  was  determined  analytically.  It  was  found  that  the 
current  is  a  consequence  of  the  formation  and  propagation  of 
interfacial  waves  at  the  thermocline.  For  large  distances  from 
the  source  of  disturbances, the  interfacial  waves  form  a  regular 
train  of  outgoing  waves  having  a  definite  amplitude  and  wave 
number.  The  surface  current  also  possesses  these  characteristic 

The  asymptotic  wave  number  of  the  surface  current  is  the 
wave  number  of  free  waves  at  the  thermocline.  It  depends  sole¬ 
ly  on  the  density  profile  and  the  frequency  of  the  disturbance 
but  not  at  all  on  its  nature  or  depth  of  submergence.  In  Fig¬ 
ure  2  the  asymptotic  wave  number,  denoted  by  a  ,  is  plotted 
against  the  frequency  uu  for  various  density  profiles.  For  a 
given  frequency  a)  and  the  same  H,  a  thicker  thermocline  retards 
the  phase  velocity,  hence  Increases  the  wave  number.  A  strati¬ 
fied  layer  below  the  thermocline,  however,  increases  the  phase 
velocity,  hence  reduces  the  wave  number.  When  the  stratifica¬ 
tion  in  the  lower  layer  is  so  large  that  its  Vaisala  frequency 

N  =V2gba  is  greater  than  the  frequency  of  the  disturbance,  inter 
facial  waves  at  the  thermocline  will  attenuate  as  they  propagate 
because  of  energy  transfer  in  the  form  of  true  internal  waves 
to  the  lower  layer. 
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The  asymptotic  amplitude  of  the  surface  current  depends  on 
how  effective  the  submerged  disturbance  is  transmitting  energy 
to  that  component  of  the  interfacial  wave  which  has  a  wave  num¬ 
ber  equal  to  a (the  wave  number  of  free  interfacial  waves 
corresponding  to  the  frequency  uu  of  the  submerged  disturbance). 
In  Figures  4  and  5  are  shown  the  asymptotic  amplitude  of  the 
surface  current  as  a  function  of  tw  for  various  density  profiles. 
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APPENDIX 

First  consider  the  surface  currents  due  to  an  oscillating 
source  of  amplitude  A  and  simple  harmonic  time  dependence  e_^U)*' 
located  at  (x  =  0,  y  =  0)  in  an  ocean  with  the  density  profile 
shown  in  Figure  lb . 

Express  the  physical  vertical  velocity  component  v(x,y,t) 
as 

1  1  -iuit 

V  =  p2V  -  p2  V  e  [ A-l] 


The  differential  equations  and  boundary  conditions  govern- 


lng  cp,  the  Fourier  transform  of  V  with  respect  to  x,  are: 

**  -k»  +[i 
ay2  l  l 

_  2gM 

UI3  1 

a3  cpx  =  0,  f  <  y  <  H 

[A -2] 

-k*  +|l  - 

dy2  t  1 

U)3  1 

cp3  -  27rAp2[6(y)  -  bs6(y)] 

for  -  "  <  y  <  f 

Ca-33 

<?1 

=  0 

at  y  =  f  +  H 

[A-4] 

Cpl 

=  cpa 

at  y  =  f 

[A-5] 
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blCpl  +  at  y  =  f 


[A -6] 


tp5  "  C  as  y  -  - 00  [A -7] 

Solving  these  equations  simultaneously  it  is  found  that 


cpi  =  C  sinh  j3i  [(f  +  H )  -  y] 


[A -8] 


where 


C  = 


l  ~Pa  f 
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Therefore 
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from  which  the  surface  current  uq  due  to  the  oscillating  source 
is  found  to  be 


+» 
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[A-13] 


The  surface  current  due  to  an  oscillating  quadrupole  of 
amplitude  M  can  easily  be  found  from  that  due  to  the  source.  It 
is 


u 
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+® 
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e _ M  I  (ba  +  Pa  )  e  a  e  da 
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[A-14] 


It  has  been  assumed  that 


e  =  P  ~  P(y  ~  f  +  H)  <  <  1  [A-15] 

P 


Let  it  be  further  assumed  that 


a  >  >  bi  or  b3 


[ A-l6] 
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so  that 


h  -  m  yr  -  ^ 
&»  -  i“i  yr-  ^ 


and  expression  [A-l4]  simplifies  to 
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[A-17] 


For  large  x,  contributions  to  the  integral  in  Equation  [A - 17 3 
come  from  the  real  poles  of  the  integrand  only.  Therefore,  let 
us  examine  the  roots  of  the  equation 


sinn 


or 
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Assume  2gbx  >  2gb2 . 

Case  I:  <  2gb3  <  2gbi 

Equation  [A-l8]  becomes 


which  has  no  real  roots.  Therefore 

uq  =  0  for  U)3  <  2gb3  <  2gbx 


Case  II:  u)a  >  2gbx  >  2gb3 

Equation  [A-l8]  becomes 

tanh  ~\Jl  -  |  a  |  H 


which  also  has  no  real  roots.  Therefore 


u  = 
o 


0 


for  itf5  >  2gbx  >  2gba 
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Case  III:  2gb2  <  u)3  <  2gbi 


Equation  [A-l8]  becomes 


(ML  _ 

ur 


which  has  an  infinite  number  of  roots 
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Introduce  dimensionless  variables 


f'  =  f/H  ,  x'  =  x/H 

U)'  =  0)/  ■\/2gb1 

‘  Na  '  =  -\f2gbl/V2gbI 


t»  =  tV2gbi 


Then 
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Write  u  as 
o 


u 

o 


2ttM 


(Ui  +  u2  +  .  . . .  ) 


Let 


lui  |  =  Ux  * 


t Us  I  =  us*  etc. 


In  Figure  6,  ux*  and  u2*  are  plotted  as  a  function  of  a» '  for 
the  special  case  b2  =  0  and  two  different  depths  of  submergence. 

From  the  calculation, it  is  found  that  the  asymptotic  sur¬ 
face  current  is  non-zero  only  when  2gb2  <  id2  <  2gbx  (assuming 
that  2gbx  >  2gb2  ) . 

When  u)a  lies  within  the  above  range,  it  is  found  that  the 
surface  current  is  in  the  form  of  a  Fourier  Series.  This  is  in 
contrast  to  calculations  for  profiles  which  possess  a  well  de¬ 
fined  thermocline.  In  the  latter  case, there  exists  only  one 
Fourier  component. 
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FIGURE  I  -  SCHEMATIC  DIAGRAM  OF  EQUILIBRIUM  DENSITY  PROFILES 


HYDRONAUTICS,  INCORPORATED 


FIGURE  3  “  CONFIGURATION  OF  QUADRUPOLE 
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FIGURE  4  -  EFFECT  OF  THERMOCLINE  THICKNESS 
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FIGURE  5-EFFECT  OF  DENSITY  GRADIENT  IN  LOWER  LAYER 
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Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Washington  7>  D.  C.  1 

2 

1  Director 

1  U.S.  Naval  Research  Laboratory 


1  Washington,  D.  C.  20360 
1  Attn:  Codes  2000  1 

1  5500  1 

1  6100  1 

1  2020  1 

1  2027  1 

5300  1 

7300  1 

4000  1 

5000  1 

7350  1 

1 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington,  D.C.  20360 
Attn:  Codes  305  1 

335  1 

1  370  1 

374  1 

421  1 

Chief,  Bureau  of  Naval  Weapons 

Department  of  the  Navy 

1  Washington,  D.C.  20360 

Attn:  Codes  RU-2  1 

RUDC-43  1 

RRRE-4  1 

Chief  of  Naval  Operations 
1  Department  of  the  Navy 
Washington,  D.C.  20360 
Attn:  Codes  0P-71B  1 

0P-07C  1 

0P-07T  1 

OP -7 17  1 

1 
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Commanding  Officer 

U.S.  Naval  Ordnance  Laboratory 

Corona,  California  1 

Commander 

U.S.  Naval  Ordnance  Test  Station 
Pasadena  Annex 
3202  E.  Foothill  Boulevard 
Pasadena  8,  California  1 

Superintendent 

U.S.  Naval  Postgraduate  School 

Monterrey,  California 

Attn:  Library  1 

Commanding  Officer  and  Director 
U.S.  Naval  Electronic  Laboratory 
San  Diego  52,  California  1 

Defense  Documentation  Center 
Cameron  Station,  5010  Duke  St. 
Alexandria,  Virginia  20 

Fluid  Mechanics  Section 
National  Bureau  of  Standards 
Washington,  D.  C.  20360 
Attn:  Dr.  0.  B.  Schubauer  1 

National  Research  Council 
Montreal  Road 
Ottawa  2,  Canada 

Attn:  Mr.  E.  S.  Turner  1 

Institute  of  Mathematical  Sciences 

New  York  University 

25  Waverly  Place 

New  York  3>  New  York 

Attn:  Prof.  J.  J.  Stoker  1 


New  York  University 

Dept,  of  Oceanography 

University  Heights 

New  York  53 j  New  York 

Attn:  Prof.  J.W.  Pierson,  Jr.  1 

The  Johns  Hopkins  University 

Baltimore  18,  Maryland 

Attn:  Prof.  0.  M.  Phillips  1 

California  Institute  of  Tech. 

Pasadena  4,  California 

Attn:  Prof.  T.Y.  Wu  1 

Director 

Scripps  Institution  of 
Oceanography 

Applied  Oceanography  Group 
University  of  California 
La  Jolla,  California  1 


Iowa  Institute  of  Hydraulic 
Research 

State  University  of  Iowa  — . 

Iowa  City,  Iowa 

Attn:  Prof.  H.  Rouse  1 

Harvard  University 
Cambridge  38,  Massachusetts 
Attn:  Prof.  G.  F.  Carrier  1 

Massachusetts  Institute  of 
Technology 

Cambridge  39>  Massachusetts 
Attn:  Prof.  J.  F.  Kennedy  1 

Stevens  Institute  of  Technology 
Davidson  Laboratory 
Hoboken,  New  Jersey 
Attn:  Mr.  D.  Savitsky  1 

Dr.  S.  J.  Lukaslk  1 
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-■3  . 


Director 

Woods  Hole  Oceanographic 
Institute 

Woods  Hole,  Massachusetts 
Dr.  B.  Zarwyn 

Airborne  Instruments  Laboratory 
Division  of  Cutler-Hammer  Inc. 
Walt  Whitman  Road,  Melville 
Long  Island,  New  York 


Commander 

U.  S.  Naval  Ordnance  U— a.). 

White  Oak,  Maryland 
1  Attn:  Dr.  S.  J.  Raff 

(DU  5-22Ji)  1 

Er,  E.  Montr*ell 
Institute  for  Defend  -Analysis 
1566  Connecticut  Aveifl  W. 

Washington,  D.  C,  1 


Mr.  Myron  J.  Block 
Block  Associates,  Inc. 

385  Putnam  Avenue 
Cambridge  39>  Massachusetts 


Director,  Special  Proj 
Department  of  the  Nat 
Washington,  D.  C.  203s 
Attn:  Code  SP-001 


•ects  Office 


O 


Director 

Hudson  Laboratories 
Dobbs  Ferry,  New  York 

Dr.  Savet 

American  Bosch  Arma  Corporation 

Arma  Division 

Garden  City,  New  York 

HYDRONAUTICS,  Incorporated 
Pindell  School  Road 
Howard  County 
Laure 1 ,  Maryland 
Attn:  Mr.  P.  Eisenberg,  Pres. 
Mr.  M.P.  Tulin,  V.  Pres. 

Dr.  L.  L.  Higgins  * 

Space  Technology  Lab.  Inc. 
Canoga  Park  Facility 
8433  Fallbrook  Avenue 
Canoga  Park,  California 


Dr  ,  I .  Katz 

Applied  Physics  Labcra"tory 
Johns  Hopkins  Univer£l'“*;j' 
Baltimore  18,  Maryland  1 

Dr.  Harvey  Brooks 
School  of  Applied  Schtmce 
Harvard  University 
Cambridge,  Massachusetts  1 

Arthur  D.  Little,  Ine, 

35  Acorn  Park; 

Cambridge  40.,  Massachi^stts 
^  Attn:  Trident  Librar/  1 

Dr.  R.  \i .  Stewart 

c/o  Pacific  Naval  Labirr-atory- 

Esquima  It, 

Vancouver,  British  CcJ^*^ia 
,  Canada  1 


Dr.  A.  V.  Hershey 
Computation  and  Exterior 
Ballistics  Laboratory 
U.S.  Naval  Weapons  Lab. 
Dahlgren,  Virginia 


Professor  M.  J.  LigntliU-H 
Imperial  College 
London,  England 


